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We review recent observational and theoretical results on the relationship be- 
tween circumstellar accretion disks and jets in young stellar objects. We then 
present a theoretical framework that interprets jets as accretion-powered, 
centrifugally driven winds from magnetized accretion disks. Recent progress 
in the numerical simulation of such outflows is described. We also discuss 
the structure of the underlying magnetized protostellar disks, emphasizing 
the role that large-scale, open magnetic flelds can play in angular momentum 
transport. 



I. INTRODUCTION 

Two of the most remarkable aspects of star formation are the pres- 
ence of disks and of energetic outflows already during the earliest phases 
of protostellar evolution. There is now strong evidence for an apparent 
correlation between the presence of outflows and of actively accreting 
disks, which suggests that there is a physical link between them. The 
prevalent interpretation is that the outflows are powered by accretion 
and that magnetic stresses mediate the inflow and outflow processes 
and eject some of the inflowing matter from the disk surfaces. If disks 
are threaded by open magnetic field lines then the outflows can take 
the form of centrifugally driven winds. Such highly collimated winds 
carry angular momentum and may, in principle, play an important role 
in the angular momentum budget of disks and their central protostars. 

This review concentrates on the developments in the study of out- 
flows and their relationship to circumstellar disks that have occurred 
since the publication of Protostars and Planets III. The reader may 
consult Konigl and Ruden in that volume and Pudritz et al. (1991) 
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for reviews of earlier work. In Sec. II we summarize the observational 
findings on outflows, disks, and magnetic fields, and their implications. 
Sec. Ill deals with the general theory of magnetized outflows and Sec. 

IV describes numerical simulations of disk-driven MHD winds. In Sec. 

V we consider the theory of magnetized protostellar disks. Our conclu- 
sions are presented in Sec. VI. 

II. OBSERVATIONAL BACKGROUND 
A. Bipolar Outflows and Jets 

Bipolar molecular outflows and narrow atomic jets are ubiquitous phe- 
nomena in protostars. There are now more than 200 bipolar CO sources 
known (see chapter by Richer et al.): they typically appear as compar- 
atively low- velocity 25 km s~^) and moderately coUimated (length- 
to-width ratios ~ 3 — 10) lobes, although several highly collimated CO 
outflows that exhibit high velocities (> 40 km s"^) near the flow axis 
have now been detected. The mass outflow rate exhibits a continuous 
increase with the bolometric luminosity of the driving source for iboi in 
the range ^ 1 — 10® Lq. Molecular outflows are present through much 
of the embedded phase of protostars and, in fact, appear to be most 
powerful and best collimated during the earliest (Class 0) protostellar 
evolutionary phase (Bontemps et al. 1996). 

The bipolar lobes are generally understood to represent ambient 
molecular material that has been swept up by the much faster, highly 
supersonic jets that emanate from the central star/disk system (see 
chapters by Eisloffel et al. and by Hartigan et al.). Jets associated 
with low-luminosity (Z/boi < 10^ Lq) young stellar objects (YSOs) have 
velocities in the range ^ 150 — 400 km s~^, large (> 20) Mach numbers, 
and can have opening angles as small as ~ 3 — 5° on scales of 10^ — 
10* AU. The inferred mass outflow rates are ~ 10"^° - 10"^ M© yr'^ 
A significant number of outflows has also been detected by optical ob- 
servations of intermediate-mass (2 < M^/Mq < 10) Herbig Ae/Be 
stars and other high-luminosity sources (Mundt and Ray 1994; Corco- 
ran and Ray 1997). The jet speeds and mass outflow rates in these 
YSOs are, respectively, a factor ~ 2 — 3 and ~ 10 — 100 higher than 
in low-iboi objects. The total momentum delivered by the jets, tak- 
ing into account both the density corrections implied by their partial 
ionization state and the long lifetimes indicated by the detection of 
pc-scale outflows, appears to be consistent with that measured in the 
associated CO outflows (e.g., Hartigan et al. 1994; Eisloffel and Mundt 
1997). A critical review of the physical mechanisms of coupling the jets 
and the surrounding gas is given in Cabrit et al. (1997). 

The momentum discharge deduced from the bipolar-outflow obser- 
vations is typically a factor ~ 10^ higher than the radiation-pressure 
thrust iboi/c produced by the central source (e.g., Lada 1985), which 
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rules out radiative acceleration of the jets. Since the bolometric lu- 
minosity of protostars is by and large due to accretion, and since the 
ratio of jet kinetic luminosity to thrust is of order the outflow speed 

10"'^ c), it follows that the jet kinetic luminosity is on average a 
fraction ~ 0.1 of the rate at which gravitational energy is liberated by 
accretion. This high ejection efficiency is most naturally understood if 
the jets are driven magnetically. 

Magnetic fields have also been implicated in jet coUimation. A par- 
ticulary instructive case is provided by HST observations of the proto- 
typical disk/jet system HH 30 (Burrows et al. 1996). The jet in this 
source can be traced to within ^ 30 AU from the star and appears as a 
cone with an opening angle of 3° between 70 and 700 AU. The narrow- 
ness of the jet indicates some form of intrinsic coUimation since external 
density gradients would not act effectively on these small scales. Mag- 
netic coUimation is a likely candidate, made even more plausible by 
the fact that the jet appears to recoUimate: its apparent opening an- 
gle decreases to 1.9° between 350 and 10^ AU. Similar indications of 
recollimation have also been found in other jets. Given that any iner- 
tial confinement would be expected to diminish with distance from the 
source, this points to the likely role of intrinsic magnetic coUimation. 

B. Connection with Accretion Disks 

The evidence for disks around YSOs and for their link to outflows 
has been strengthened by a variety of recent observations. These in- 
clude systematic studies of the frequency of disks by means of infrared 
and millimeter surveys and further interferometric mappings (now com- 
prising also the submillimeter range) that have resolved the structure 
and velocity field of disks down to scales of a few tens of AU (see chap- 
ter by Wilner and Lay). High-resolution images of disks in several jet 
sources have also been obtained in the near-infrared using adaptive op- 
tics and in the optical using the HST (see chapter by McCaughrean et 
al.). 

Stellar jets are believed to be powered by the gravitational energy 
liberated in the accretion process and to be fed by disk material. This 
picture is supported by the strong apparent correlation that is found 
(e.g., Cabrit et al. 1990; Cabrit and Andre 1991; Hartigan et al. 1995) 
between the presence of outflow signatures (such as P Cyg line pro- 
files, forbidden line emission, thermal radio radiation, or well-developed 
molecular lobes) and accretion diagnostics (such as ultraviolet, infrared, 
and millimeter-wavelength emission excesses, or inverse P Cyg line pro- 
files). Further support is provided by the apparent decline in outflow 
activity with stellar age, which follows the similar trend exhibited by 
the disk frequency (see chapters by Andre et al. and by Mundy et 
al.) and mass accretion rate (see chapter by Calvet et al.). While 
virtually every Class source has an associated outflow, a survey of 
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optical and molecular outflows in the Taurus- Auriga cloud (Gomez et 
al. 1997) found an incidence rate of ^ 60% among Class I objects but 
only ~ 10% among Class II ones (and none in Class III objects). The 
inference that jets are powered by accretion and originate in accretion 
disks is strengthened by the evidence for disks in the youngest YSOs 
in which outflows are detected. For example, sub-mm interferometric 
observations of VLA 1623, one of the youngest known Class sources, 
imply the presence of a circumstellar disk of radius < 175 AU and mass 
> .03 M0 (Pudritz et al. 1996). 

Corcoran and Ray (1998) demonstrated that the correlation be- 
tween [OIJA6300 line luminosity (an outflow signature) and excess in- 
frared luminosity (an accretion diagnostic) originally found in Class 
II sources extends smoothly to YSOs with masses of up to ^ 10 Mq 
and spans 5 orders of magnitude in luminosity. It is noteworthy that 
correlations of the type M oc L^qj that apply to both low-luminosity 
and high-luminosity YSOs have been established in several independent 
studies for the mass accretion rate (from IR continuum measurements; 
Hillenbrand et al. 1992 - see, however, Hartmann et al. 1993, Bell 
1994, Miroshnichcnko et al. 1997, and Pezzuto et al. 1997 for alter- 
native interpretations of the infrared emission in Herbig Ae/Be stars), 
the ionized mass outflow rate in the jets (from radio continuum obser- 
vations; Skinner et al. 1993), and the bipolar molecular outflow rate 
(from CO line measurements; Levreault 1988). Taken together, these 
relationships suggest that a strong link between accretion and outflow 
exists also in high-mass YSOs and that the underlying physical mech- 
anism is basically the same as in low- mass objects (see Konigl 1999). 

Strong evidence for a disk origin of jets is available for the energetic 
outflows associated with FU Orionis outbursts (see chapters by Bell et 
al. and by Calvet et al.). The outbursts have been inferred to arise in 
young YSOs that are still rapidly accreting, although it is possible that 
they last into the Class II phase. The duration of a typical outburst 
is ~ 10^ yr, and during that time the mass accretion rate (as inferred 
from the bolometric luminosity) is ~ 10~^ Mq yr~^, with the deduced 
mass outflow rate A^wind (at least in the most powerful sources like FU 
Ori and Z CMa) being a tenth as large. The ratio Mwind/Afacc ~ 0.1 is 
similar to that inferred in Class II YSOs and again points to a rather 
efficient outflow mechanism. Detailed spectral modeling demonstrates 
that virtually all the emission during an outburst is produced in a ro- 
tating disk. Furthermore, the correlation found in the prototype FU 
Ori between the strength and the velocity shift of various photospheric 
absorption lines can be naturally interpreted in terms of a wind accel- 
erating from the disk surface (Calvet et al. 1993; Hartmann and Calvet 
1995). Because of the comparatively low temperatures (^ 6000 K) in 
the wind acceleration zone, thermal-pressure and radiative driving are 
unimportant (the latter due to the lack of a high-temperature source): 
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magnetic driving is thus strongly indicated. The recurrence time of 
outbursts has been estimated to he in the range ~ 10"^ — 10^ yr, and if 
these outbursts are associated with the large-scale bow shocks detected 
in pc-scale jets (e.g., Reipurth 1991), then a value near the lower end 
of the range is implied. In that case most of the stellar mass would be 
accumulated through this process, and, correspondingly, most of the 
mass and momentum ejected over the lifetime of the YSO would orig- 
inate in a disk-driven outflow during the outburst phases (Hartmann 
1997). 

C. Magnetic Fields in Outflow Sources 

The commonly accepted scenario for the origin of low-mass proto- 
stars is that they are produced from the collapse of the inner regions of 
molecular clouds that are supported by large-scale magnetic fields (and 
likely also hydromagnetic waves). In this picture, a gravitationally un- 
stable inner core forms as a result of mass redistribution by ambipolar 
diffusion and subsequently collapses dynamically (see review by McKee 
et al. 1993). The mass accretion rates predicted by this picture are 
consistent with the inferred evolution of yoimg YSOs (e.g., Ciolck and 
Konigl 1998). Basic support for this scenario is provided by far- infrared 
(e.g., Hildebrand et al. 1995) and sub- mm (e.g.. Greaves et al. 1994, 
1995; Schleiming 1998; Greaves and Holland 1998) polarization mea- 
surements, that reveal an ordered, hourglass-shaped field morphology 
on sub-pc scales consistent with the field lines being pulled in at the 
equatorial plane of the contracting core. Moreover, H I and OH Zceman 
measurements (e.g., Crutcher et al. 1993, 1994, 1996) are consistent 
with the magnetic field having the strength to support the bulk of the 
cloud against gravitational collapse. 

There now exist measurements of magnetic fields in the flows them- 
selves at large distances from the origin (see chapter by Eisloffel et al.). 
In particular, the strong circular polarization detected in T Tau S in 
two oppositely directed nonthermal emission knots separated by 20 AU 
indicates a field strength of at least several gauss (Ray et al. 1997). 
This high value can be attributed to a magnetic field that is advected 
from the origin by the associated stellar outflow and that dominates the 
interal energy of the jet. This observation thus provides direct evidence 
for the essentially hydromagnetic character of jets. 

III. MHD WINDS FROM ACCRETION DISKS 
A. Basic MHD Wind Theory 

The theory of centrifugally driven winds was first formulated in the 

context of rotating, magnetized stars (Schatzmann 1962; Weber and 
Davis 1967; Mestel 1968). Using 1-D, axisymmetric models, it was 
shown that such stars could lose angular momentum by driving winds 
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of this type. This idea was apphcd to magnetized accretion disks in 
the seminal paper of Blandford and Payne (1982, BP). Every annulus 
of a Keplerian disk may be regarded as rotating close to its "breakup" 
speed, so disks arc ideal drivers of outflow when sufficiently well mag- 
netized. The removal of disk angular momentum allows matter to move 
inward and produces an accretion flow. In a steady state, field lines 
must also slip radially out of the accreting gas and maintain their fixed 
position in space: this necessitates diffusive processes. As we discuss in 
Sec. V.B, strong field diffusivity is a natural attribute of the partially 
ionized regions of protostellar disks, and it could counter both the ad- 
vection of the field lines by the radial inflow and their winding-up by 
the differential rotation in the disk. 

Consider, for the moment, the simplest possible description of a 
magnetized, rotating gas threaded by a large-scale, open field (charac- 
terized by an even symmetry about the midplane z = 0). The equations 
of stationary, axisymmetric, ideal MHD arc the conservation of mass 
(continuity equation); the equation of motion with conducting gas of 
density p subject to forces associated with the pressure P, the gravi- 
tational field (from the central object whose gravitational potential is 
(/>), and the magnetic field B; the induction equation for the evolution 
of the magnetic field in the moving fluid; and the solenoidal condition 
on B: 

V-(pV) = 0, (1) 

pV-VV = -Vp-pV(/)+-^(VxB)xB, (2) 

47r 

V X (V X B) = , (3) 



V • B = . (4) 

Consider the angular momentum equation for axisymmetric 
flows. This is described by the cp component of equation (2). Ignoring 
stresses that would arise due to turbulence, and noting that neither the 
pressure nor the gravitational term contributes, one finds that 

pVp • V{rV^) = ^ • V(rB^) , (5) 

where we have broken the magnetic and velocity fields into poloidal 
and toroidal components: B = Bp -h S^e^ and V = Vp -I- V^e^. 

Important links between the velocity field and the magnetic field 
are contained in the induction equation (3), whose solution is 



V X B = W, 



(6) 
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where ip is some scalar potential. This shows that the electric field due 
to the bulk motion of conducting gas in the magnetic field is derivable 
from an electrostatic potential. This has two important ramifications. 
The first is that, because of axisymmetry, the toroidal component of 
this equation must vanish {dip/d4> = 0). This forces the poloidal ve- 
locity vector to be parallel to the poloidal component of the magnetic 
field, Vpl |Bp. This, in turn, implies that there is a function k, the mass 
load of the wind, such that 

pVp = /cBp . (7) 

Substitution of this result into the continuity equation (1), and then 
use of the solenoidal condition (eq. [4]) reveals that is a constant 
along a surface of constant magnetic flux, i.e., that it is conserved 
along field lines. This function can be more rcvcalingly cast by noting 
that the wind mass loss rate passing through an annular section of the 
flow of area dA through the flow is dM^ = pVpdA, while the amount 
of poloidal magnetic flux through this same annulus is d^ = BpdA. 
Thus, the mass load per unit time and per unit magnetic flux, which 
is preserved along each streamline emanating from the rotor (a disk in 
this case), is 

pFp _ dM^ 
"-^-^^ 

The mass load is determined by the physics of the underlying rotor, 
which is its source. 

A second major consequence of the induction equation follows from 
the poloidal part of equation (6). Taking the dot product of it with 
Bp and using equation (7), one easily proves that ij) is also a con- 
stant along a magnetic fiux surface and that it must take the form 
Ip = fl — {kB^/pr). In order to evaluate V, note that = at the 
disk midplane by the assumed even symmetry. Thus t]j equals fio, the 
angular velocity of the disk at the midplane. One thus has a relation 
between the toroidal field in a rotating flow and the rotation of that 
flow, 

B4, = ^{n-no) . (9) 

Let us now examine the angular momentum equation. Returning 
to the full equation (5) and applying equation (7) and the constancy 
of k along a field line, one obtains 

Bp • V(ry, - ^) = . (10) 
Hence the angular momentum per unit mass, 
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is constant along a streamline. This shows that the specific angular 
momentum of a magnetized flow is carried by both the rotating gas 
(first term) and the twisted field (second term). The value of I may be 
found by eliminating the toroidal field between equations (9) and (11) 
and solving for the rotation speed of the flow, 

; 2 _ 2q2 

rV, = ^HL^Llk , (12) 
m"^ — 1 

where the Alfvcn Mach number m of the flow is defined as = /V^, 
with Va. = Bp/ [A.-npf-l'^ being the Alfven speed of the flow. The Alfven 
surface is the locus of the points r = ta on the outflow field lines where 
m = 1. The fiow along any field line essentially corotates with the 
rotor until this point is reached. ^From the regularity condition at the 
Alfven critical point (where the denominator of eq. [12] vanishes) one 
infers that the conserved specific angular momentum satisfies 

I = ^orl . (13) 

If we imagine following a field line from its footpoint at a radius vq, 
the Alfven radius is at a distance rA(ro) from the rotation axis and 
constitutes the lever arm for the back torque that this flow exerts on 
the disk. The other critical points of the outflow are where the outflow 
speed Vp equals the speed of the slow and fast magnetosonic modes in 
the flow (at the so-called SM and FM surfaces). 

Finally, a generalized version of Bernoulli's equation may be de- 
rived by taking the dot product of the equation of motion with Bp. 
One then finds that the specific energy 

E = ^(Vp + ^^r^) +<l> + h + f^o(J^ori - flr^) , (14) 

where h is the enthalpy per unit mass, is also a field- line constant. 

The terminal speed Vp = Voc corresponds to the region where the 
gravitational potential and the rotational energy of the flow are neg- 
ligible. Since for cold fiows the specific enthalpy may be ignored, one 
infers from equation (14) 

~ , (15) 

a result first obtained by Michel (1969) for 1-D fiows. The impor- 
tant point regarding outflow speeds from disks is that K»/i^o^o ~ 
2^/^rA/7'o: the asymptotic speed is larger than the rotor speed by a 
factor that is approximately the ratio of the lever arm to the footpoint 
radius. 
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B. Connection with Underlying Accretion Disk 

We now apply the angular momentum conservation relation (5) to cal- 
culate the torque exerted on a thin accretion disk by the external mag- 
netic field. The vertical flow speed in the disk is negligible, so only the 
radial inflow speed Vr and the rotation speed (Keplerian for thin 
disks) contribute. On the right-hand side, both the radial and vertical 
magnetic contributions come into play, so 

ro dro Atztq dro 47r dz 

One sees that specific angular momentum is removed from the inward 
accretion flow by the action of two types of magnetic torque. The 
first term on the r.h.s. represents radial angular momentum associated 
with the radial shear of the toroidal field, while the second term is 
vertical tranport due to the vertical shear of the toroidal field. In a thin 
disk, and for typical field inclinations, the second term will dominate. 
Note that the first term vanishes at the disk midplanc because i?r = 
there. Now, following standard thin-disk theory, vertical integration of 
the resulting equation gives a relation between the disk accretion rate, 
-^acc = —2TT'EVrro, and the magnetic torques acting on its surfaces 
(subscript s), 

JWacc -J^ — — rg±!0,si>z , [i-l) 

Angular momentum is thus extracted out of disks threaded by open 
magnetic fields. The angular momentum can be carried away either by 
toroisional Alfvcn waves or, when the magnetic field lines are inclined 
by more than 30° from the vertical, by a ccntrifugally driven wind. By 
rewriting equation (11) as ri?^ = 47rfc(rV0 — I) and using the derived 
relations for k and I, the disk angular momentum equation can be cast 
into its most fundamental form, 

^^JJMI = 'J^n,rlil - ir,/r,r) . (18) 
ara dro 

This equation shows that there is a crucial link between the mass out- 
fiow in the wind and the mass accretion rate through the disk: 

Mace ^ {rA/rof Mwind • (19) 

One has arrived at the profoundly useful expression of the idea that, 
if viscous torques in the disk are relatively unimportant, the rate at 
which the disk loses angular momentum {jd = Macc^o^o) exactly 
the rate at which it is carried away by the wind [j^ = M^ind^o'i^A) ■ 

The value of the ratio rx/ro is ~ 3 for typical parameters, so 
one finds Mwind/^facc — 0.1, which is in excellent agreement with the 
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observations (Sec. II). The explanation of this relationship is thus 
intimately linked to the disk's angular momentum loss to the wind. 

C. Flow Initiation and Collimation 

Stellar winds usually require hot coronae to get started, while winds 
from disks, which effectively rotate near "breakup," do not. In the 
case of a thin disk even a cool atmosphere will suffice as long as the 
field lines emerging from the disk make an angle < 60° to the surface. 
This follows from Bernoulli's equation by comparing the variations in 
the effective gravitational potential and the kinetic energy of a particle 
that moves along a field line near the disk surface (BP; Konigl and 
Ruden 1993; Spruit 1996; but see Ogilvie and Livio 1998). 

The collimation of an outflow, as it accelerates away from the disk, 
arises to a large extent from the hoop stress of the toroidal field compo- 
nent. Prom equation (9) one sees that at the Alfven surface jB^j ~ Bp 
and that in the far field (assuming that the flow opens up to radii 
r » ta) -B^/Bp — r/rA- Thus, the inertia of the gas, forced to 
corotate with the outflow out to ta, causes the jet to eventually self- 
collimate through the jz x B0 force, which is known as the z-pinch in 
the plasma physics literature. 

The detailed radial structure of the outflow is deduced by balanc- 
ing all forces perpendicular to the field lines and is described by the 
so-called Grad-Shafranov equation. This is a complicated nonlinear 
equation for which no general solutions are available. Because of the 
mathematical difficulties (e.g., Hcinemann and Olbert 1978), the ana- 
lytic studies have been characterized by simplified approaches, includ- 
ing separation of variables (e.g., Tsinganos and Trussoni 1990; Sauty 
and Tsinganos 1994), self-similarity (e.g., BP; Bacciotti and Chiuderi 
1992; Contopoulos and Lovelace 1994; Lynden-BeU and Boily 1994), 
previously "guessed" magnetic configurations (e.g., Pudritz and Nor- 
man 1983; Lery et al. 1998), as well as examinations of various asymp- 
totic limits to the theory (e.g., Heyvaerts and Norman 1989, HN; Appl 
and Camenzind 1993; Ostriker 1997). Pelletier and Pudritz (1992) con- 
structed non self-similar models of disk winds including cases where the 
wind emerges only from a finite portion of the disk. 

Self-similarity imposes a specific structure on the underlying disk. 
In the BP model, all quantities scale as power laws of spherical radius 
along a given direction. This directly implies that the Alfven surface is 
conical and, similarly, that the disk's scale height H{r) scales linearly 
with r. Furthermore, inasmuch as the problem contains a characteristic 
speed, namely the Kepler speed in the disk, one infers that Va cx Cg oc 
V[ oc oc Vk, i.e., that the Alfven, sound, radial inflow, terminal 
outflow, and Kepler speeds, respectively, arc all proportional to one 
another. The scaling Cg oc Vk implies that the disk temperature has the 
virial scaling T oc r~^. The scaling 14 a Vk as well as the H{r) relation 
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imply that, for a constant mass accretion rate Afacc = —2T:{2Hp)Vrr, 
the density p oc r~'^/^. Next, the scahng of the disk Alfven speed 
Va oc Vk together with the density result imply that the disk poloidal 
field (and hence also B^) scales as r~^^^. In turn, the mass load fco oc 

Tq^^'' and the wind mass loss rate Mwind(?'o) Inro. More general 
self-similar models may be constructed by making the more realistic 
assumption that mass is lost from the disk so that the accretion rate is 
not constant; Macci''') oc r~'^. Adopting the radial scaling of the disk 
magnetic field, Bo{r) oc r~'^, one finds that self-similarity imposes the 
scaling /x = 2(z^ — 1.25) (see eq. [17]). 

The work of HN and others (see Heyvaerts and Norman, 1997) has 
shown that in general two types of collimation are possible, depending 
on the asymptotic behavior of the current intensity / = {c/2)rB^ oc 
r^pfio- If I ^ as r — > oo, then the field lines are space-filling 
paraboloids, whereas if this limit for the current is finite, then the flow 
is collimated to cylinders. The character of the flow therefore depends 
on the boundary conditions at the disk. 

Finally, we note that the stability of magnetized jets with toroidal 
magnetic fields is strongly assisted by the jet's poloidal magnetic field, 
which acts as a spinal column for the jet (e.g., Appl and Camenzind 
1992). 

IV. NUMERICAL SIMULATIONS OF DISK WINDS 

The advent of numerical simulations has finally made it possible 
to study the rich, time-dependent behavior of MHD disk winds. This 
allows one to test the stationary theory presented above as well as to 
search for the conditions that give rise to episodic outflows. These 
simulations represent perhaps the main advance in the subject since 
Protostars and Planets III. The published simulations generally assume 
ideal MHD and may be grouped into two classes: 1) dynamic MHD 
disks, in which the structure and evolution of the magnetized disk is 
also part of the simulation, and 2) stationary MHD disks, in which the 
underlying accretion disk does not change and provides fixed boundary 
conditions for the outflow problem. 

A. Dynamic Disks and Winds 

The flrst numerical calculations of disk winds were published by Uchida 
and Shibata (1985, US), and Shibata and Uchida (1986). These simula- 
tions modeled a magnetized disk in sub-Keplerian rotation and showed 
that a rapid radial collapse develops in which the initially poloidal field 
threading the disk is wound up due to the differential rotation. The ver- 
tical Alfven speed, being smaller than the free-fall speed, implies that 
a strong, vertical toroidal field pressure gradient dB'^/dz must rapidly 
build up. This force results in the transient ejection of coronal material 
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above and below the disk as the spring uncoils. The work of US was 
confirmed by Stone and Norman (1994, SN) using their ideal-MHD, 
ZEUS 2-D code (Stone and Norman 1992). 

ff the mean magnetic field energy density is less than the thermal 
pressure in the disk {Va < Cg), then a strong magnetorotational in- 
stability will develop (e.g., Balbus and Hawley 1991, BH). In 2-D, this 
leads to a vigorous, radial, channel flow and rapid outward transport of 
angular momentum. SN ran a series of ZEUS 2-D simulations for a uni- 
form magnetic field threading a wedge-shaped disk and a surrounding 
corona (models defined by 4 parameters). Their simulations investi- 
gated three cases (see also Bell and Lucek 1995 and Matsumoto et al. 
1996): (A) sub-Keplerian rotation, (B) a Kepler disk with strong disk 
field, and (C) a Kepler disk with a weak field. In Case A, rapid collapse 
immediately ensues with an expanding, transient outflow appearing in 
2.5 orbits (reproducing US). In Case B, the disk is BH-stable, but col- 
lapse occurs anyway because of the very strong braking of the disk due 
to the external MHD torque. In case C one again sees rapid radial 
collapse of the disk, this time because of a strong BH instability. 

The 2-D channel flow does not, however, persist in 3-D: in that case 
the BH instability develops into a fully turbulent flow and the inflow 
rate is significantly reduced (see chapter by Stone et al.). 

B. Stationary Disks and Winds 

The outflow problem can be clarified by focusing on the more restricted 
question of how a wind is accelerated and coUimated for a prescribed 
set of fixed boundary conditions on the disk. Keplerian disks in 3-D are 
stable on many tens of orbital times and this justifies the stationary disk 
approach: the launch and coUimation of jets from their surfaces occurs 
in only a few inner-disk rotation times. Groups that have taken this 
route include Ustyugova et al. (1995), Ouyed et al. (1997, OPS), Ouyed 
and Pudritz (1997a,b, OPI and CPU), Romanova et al. (1997), and 
Meier et al. (1997). The published simulations differ in their assumed 
initial conditions, such as the magnetic field distribution on the disk, 
the plasma /? (= Pgas/-Pmag) above the disk surfaces, the state of the 
initial disk corona, and the handling of the gravity of the central star. 
Broadly speaking, all of the existing calculations show that winds from 
accretion disks can indeed be launched and accelerated, much along the 
lines suggested by the theory presented in Sec. HI. The results differ, 
however, in the degree to which flow collimation occurs. 

Ustyugova et al. (1995) and Romanova et al. (1997) employed 
a magnetic configuration described by a monopole field centered be- 
neath the disk surface. Nonequilibrium initial conditions as well as a 
softened gravitational potential were used. Relatively low resolution 
simulations of flows with /? >> 1 (Ustyugova et al.) showed that col- 
limated, non-stationary outflows develop. Similar simulations in the 
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strong-field (/? < 1) regime (Romanova et al.) resulted in stationary, 
but uncoUimated outflows on these scales. 

Simulations by OPS, OPI, and OPII (see Pudritz and Ouyed 1997 
for a review) employed the ZEUS 2-D code and studied two different 
magnetic configurations. The coronal gas was initially taken to be in 
hydrostatic equilibrium with the central object as well as in pressure 
balance with the top of the disk. Two initial magnetic configurations 
were adopted, each chosen so that no Lorentz force is exerted on the 
hydrostatic corona, i.e., such that j = 0. These are the potential field 
configuration of Cao and Spruit (1994), and a uniform, vertical field 
that is everywhere parallel to the disk rotation axis. The gravity was 
unsoftened in these calculations, (500 x 200) spatial zones were used, 
and simultations were run up to 400ti (where tj is the Kepler time for 
an orbit at the inner edge of the disk, r^). This model is described by 
5 parameters set at ri: 3 to describe the initial corona (e.g., (3, = 1.0), 
as well as 2 parameters to describe the disk physics (e.g., the injection 
speed Vinj of the material from the disk into the base of the corona). 
The initial conditions correspond to turning on the rotation of the 
underlying Keplerian disk at t = 0. 

The first thing that happens is the launch and propagation of a brief 
transient, torsional Alfven wave front that sweeps out from the disk 
surface but that leaves the corona largely undisturbed. An outflow also 
begins almost immediately and develops into a stationary or an episodic 
jet. These jet-like outflows are highly collimated and terminate in a jet 
shock. A bow shock, driven by the jet, pushes through the corona. A 
noticeably empty cavity dominates most of the volume behind the bow 
shock: the cavity is filled by the toroidal magnetic field that is generated 
by the jet itself. This may provide an explanation of the extended 
bipolar cavities that surround highly collimated jets. For a fiducial 
injection speed of 10~^Vk, the potential field configuration developed 
into a stationary outflow with many (but not all, since the flow is 
not self-similar) of the characteristics of the BP solution. The Alfven 
surface is correctly predicted by the analysis in section III. The Alfven 
and fast-magnetosonic Mach numbers reach 5 and 1.6, respectively, at 
z = lOn, with the toroidal-to-poloidal field-strength ratio being ~ 3 on 
this scale. Outflow only takes place on field lines that are inclined by 
less than 60° with respect to the disk surface, as predicted by BP. 

Figure 1 (adapted from OPS) shows that the poloidal field lines in 
the initial state are collimated toward the rotation axis by the hoop 
stress of the jet's toroidal field. The figure also shows the propagation 
of the jet-driven bow-shock and the eventual creation of a cylindri- 
cally collimated jet with well determined Alfven and fast-magnetosonic 
critical surfaces in the acceleration region above the disk (the slow- 
magnetosonic surface is too close to the disk to be resolved in this 
figure). 
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Figure 1. Numerical sirinilations of disk-driven MHD outflows (adapted from 
OPS). In frame a, the left panel shows the initial magnetic configuration cor- 
responding to the "potential field" solution, whereas the right panel displays 
the initial isodensity contours of the corona. The flow injection speed is 
Vinj = 10~^Vk. Frames h and c show the evolution of the initial magnetic 
and density structures (the left and right panels, respectively) at 100 and 400 
inner time units. Frame c also displays the locations of the Alfvcn critical 
surface (filled hexagons) and the FM surface (stars). In frame d, the Alfven 
lever arm (rA/ro) found in the simulation is shown (filled hexagons) as a 
function of the location (ro) of the footpoints of the field lines on the disk, 
and is compared to the prediction from the steady-state theory of Sec. III. A 
(squares) . 



For the same boundary and initial conditions, the initially vertical 
field configuration leads to the development of a jet that is episodic 
over the 400 ti duration of the simulation (OPS, OPII). Even though 
the initial configuration is highly unfavorable to jet formation, jet pro- 
duction occurs. The reason for this is the effect of the toroidal field in 
the corona, which is concentrated towards the inner edge of the disk 
(where the Kepler rotation is the greatest). It therefore exerts a radial 
pressure force d{B-^/8n)/dr that opens up field lines at larger radii. 
As long as /3 ~ 1 in this region, field lines are pliable enough to move, 
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and a jet is launched. Episodic knots are produced on a time; sc;ale 
Tknot — fjet/VA,<i>, which is reminiscent of a type of kink instability. 
Regions of high toroidal field strength and low density separate and 
confine the knots, which in turn have high density and low toroidal 
field strength. 

The appearance of stationary or episodic jets has nothing to do 

with the magnetic configuration, but rather with the mass loading of 
the magnetic configuration, as is seen in an extensive set of simulations 
(Ouyed and Pudritz 1998). Figure 2 shows a simulation with the same 
potential configuration parameters as in Figure 1 except for a reduced 
injection speed (and hence mass load) of Vinj = 10~^Vk- Clear episodic 
behavior is now seen in this flow. 
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Figure 2. MHD outflow simulations with the same paramotors and magnetic 
configuration as in Fig. 1, except Vinj = 10~^Fk (adapted from Ouyed and 
Pudritz 1998). The left panels show the magnetic field structure and knot 
formation at times 0.0, 75.0, 112.5, and 150.0 ti. The right panels show the 
angle 9o between the field lines and the disk surface at these 4 times. Only 
field lines that open up to an angle 6o < 60° are seen to drive an outflow. 



Nonstcady outflows could conceivably also arise in highly conduct- 
ing disk regions where the winding-up of the field lines might eject gas 
through a strong \B^\ pressure gradient (e.g., Contopoulos 1995). 
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C. Far-Field Behavior of MHD Jets 

Simulations of MHD jets propagating into a uniform medium show that, 
in certain respects, they differ dramatically from hydrodynamic jets. 
Clarke et al. (1986), Lind et al. (1986), and Kossel et al (1990) showed 
that the strong toroidal field of a jet prevents matter from spraying side- 
ways on encountering the jet shock. Rather, jet material decelerated by 
a Mach disk and a strong annular shock is focused mainly forward into 
a nose cone. This contrasts with the backflowing cocoon that charac- 
terizes purely hydrodynamic, low-density jets. Hydromagnetic jets are 
generally preceeded by stronger, more oblique shocks and advance far 
more quickly into the surrounding medium than their hydrodynamic 
counterparts (see Kossel et al. 1990 for more detail). MHD jets may 
thus explain why one does not measure strong transverse expansions in 
bipolar molecular lobes (e.g., Masson and Chernin 1993; Cabrit 1997). 

V. MAGNETIZED PROTOSTELLAR ACCRETION DISKS 

A. Magnetic Angular Momentum Transport 

One of the primary reasons why open, ordered magnetic fields are im- 
portant in accreting astrophysical systems is that they can mediate 
a vertical transport of angular momentum from accretion disks. In 
Sec. HI we discussed the angular momentum transport by centrifugally 
driven winds, but it is important to realize that angular momentum can 
be removed from the disk surfaces even if the ficld-linc inclination with 
respect to the vertical is not large enough for the wind launching con- 
dition to be satisfied, so long as the field is attached to a "load" that 
exerts a back torque on the disk. In that case the disk can lose angu- 
lar momentum through torsional Alfven waves propagating away from 
the disk surfaces (the "magnetic braking" mechanism; see Mouschovias 
1991 for a review). When the magnetic fields are well below cquiparti- 
tion with the gas pressure in a differentially rotating disk, the magne- 
torotational (Balbus-Hawley) instability develops on a dynamical time 
scale and produces a magnetic stress-dominated turbulence that gives 
rise to radial transport of angular momentum characterized by an ef- 
fective viscosity parameter a that lies in the range 0.005 — 0.5 (see 
chapter by Stone et al.). 

Numerical simulations (see Sec. IV. A) have demonstrated that 
all of the above processes could contribute to the angular momentum 
transport in magnetized disks. One can estimate the relative roles of 
vertical wind transport and radial turbulent transport by taking the 
ratio of the external wind torque computed in Sec. Ill and the viscous 
torque associated with magnetoccntrifugal turbulence, which can be 
written as {Bl/ AirP^^) {ta/ aH) (Pelletier and Pudritz 1992). Thus, 
even if the ordered field threading the disk were far below equipartition 
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with the disk gas prc;ssurc, a wind torque could stiU dominate a tur- 
bulent torque simply because of its large lever arm (the Alfven radius 
ta, which greatly exceeds the disk scale height H). We note in passing 
that a long effective lever arm also characterizes angular momentum 
transport by spiral waves in the disk. 

The requirement that, in a steady state, the torque exerted by a 
large-scale magnetic field at the surface of a thin and nearly Keplerian 
accretion disk balances the inward angular momentum advection rate 
can be written as Mace = 2r^/2 |Bz B^,s|/(G'M»)i/^ (see eq. [17]). 
Assuming a rough equality between the vertical and azimuthal surface 
field components, this relation implies that, at a distance of 1 AU from 
a solar-mass protostar, a 1 G field (the value indicated by meteoritic 
data for the protosolar nebula; Levy and Sonett 1978) could induce 
accretion at a rate of ~ 2 x 10~^ Mq yr'^, which is compatible with 
the mean values inferred in embedded protostars (see Sec. II. B). For 
comparison, the minimum-mass solar nebula model implies a thermal 
pressure at 1 AU that corresponds to an equipartition magnetic field 
of 18 G, which implies that a ~ 1 G field could be readily anchored 
in the protosolar disk at that location. 

B. Wind-Driving Protostellar Disk Models 

The properties of protostellar accretion disks and of integrated 
disk/wind systems have, so far, been derived only under highly sim- 
plified assumptions. Among the papers that can be consulted on this 
topic are Konigl (1989, 1997), Wardle and Konigl (1993), Ferreira and 
Pelletier (1993a,b, 1995; also Ferreira 1997), Lubow et al. (1994), Li 
(1995, 1996), and Reyes- Ruiz and Stepinski (1996). A plausible ori- 
gin for an open disk magnetic field is the interstellar field that had 
originally threaded the magnetically supported molecular cloud and 
that was subsequently carried in by the collapsing core. This picture 
is favored over disk-dynamo interpretations in view of the fact that 
the latter typically produce closed (quadrupolar) field configurations, 
although we note that scenarios for opening dynamo-generated field 
lines have been considered in the literature (e.g.. Tout and Pringle 
1996; Curry, Pudritz, and Sutherland 1994). 

The theory and simulations of the previous sections assumed ideal 
MHD. While this approximation is adequate for the surface layers of 
protostellar disks, the disk interiors are typically weakly ionized, and 
their study entails the application of multifluid MHD. For typical pa- 
rameters of disks around solar-mass YSOs, one can distinguish between 
the low-density regime (on scales r ^ 100 AU), in which the current 
in the disk is carried by metal ions and electrons (whose densities are 
determined from the balance between ionizations by cosmic rays and re- 
combinations on grain surfaces) , and the high-density regime (on scales 
r ^1 — 10 AU), where the current is carried by small charged grains or 
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by ions and electrons that rccombinc without grains. We restrict our 
attention to the low-density regime and concentrate on the case where 
the magnetic field B is "frozen" into the electrons and diffuses relative 
to the dominant neutral component as a result of an ion-neutral drift 
(ambipolar diffusion). 

To derive a self-consistent steady-state disk/wind configuration, 
one combines the mass, momentum (radial, vertical, and angular), and 
energy conservation relations, together with Maxwell's equations and 
the generalized Ohm's law, and imposes the requirements that the out- 
flow pass through the relevant critical points (sec Sec. IV). So far 
only simple prescriptions for the disk thermal structure (isothermal or 
adiabatic) and conductivity (ambipolar diffusion in the density regime 
where the ion density is constant or Ohmic diffusivity parametrized 
using a "turbulence" prescription) have been considered. These are 
probably adequate for obtaining the basic structure of the disk, but 
more realistic calculations are needed to correctly model the transition 
region between the disk and the wind. 

The vertical structure of a generic centrifugal wind-driving disk 
(or "active" surface layer; see Sec. V.C) can be divided into 3 distinct 
zones (see Fig. 3): a quasi- hydrostatic region near the midplane of 
the disk, where the bulk of the matter is concentrated and most of the 
field-line bending takes place, a transition zone where the inflow grad- 
ually diminishes with height, and an outflow region that corresponds 
to the base of the wind. The first two regions are characterized by a 
radial inflow and sub-Keplerian rotation, while the gas at the base of 
the wind flows out with Vrj, > Vk- The boundary conditions for the sta- 
tionary disk simulations discussed in Sec. IV. B arise from the physical 
properties of this latter region. 

What determines the disk structure and how does the magnetic 
field extract the angular momentum of the accreting gas? The quasi- 
hydrostatic region is matter dominated, with the ionized plasma and 
magnetic field being carried around by the neutral material. The ions 
are braked by a magnetic torque, which is transmitted to the neutral gas 
through the frictional (ambipolar-diffusion) drag; therefore Vi0 < in 
this region (with the subscript i denoting ions). The neutrals thus 
lose angular momentum to the field, and their back reaction leads to a 
buildup of the azimuthal field component [S^j away from the midplane. 
The loss of angular momentum enables the neutrals to drift toward the 
center, and in doing so they exert a radial drag on the field lines. This 
drag must be balanced by magnetic tension, so the field lines bend away 
from the rotation axis. This bending builds up the ratio B^/B^, which 
needs to exceed l/\/3 at the disk surface to launch a centrifugally 
driven wind. The magnetic tension force, transmitted through ion- 
neutral collisions, contributes to the radial support of the neutral gas 
and causes it to rotate at sub-Keplerian speeds. 
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neutrals 

^ >- ions 

Figure 3. Schematic diagram of the vertical structure of an ambipolar 
diffusion-dominated disk, showing a representative field line and the 
poloidal velocities of the neutral (solid arrowheads) and the ionized (open 
arrowheads) fluid components. Note that the poloidal velocity of the ions 
vanishes at the midplane {z = 0) and is small for both fluids at the base of 
the wind (z = Zb). The relationship between the azimuthal velocities is also 
indicated. 



The growth of the radial and azimuthal field components on moving 
away from the midplane results in a magnetic pressure gradient that 
tends to compress the disk. The vertical compression by the combined 
magnetic and tidal stresses is, in turn, balanced by the thermal pressure 
gradient. The magnetic energy density becomes dominant as the gas 
density decreases, marking the beginning of the transition zone (at z ~ 
Zh). The field above this point is nearly force free ([V x B]xB « 0), 
so the field lines, which vary only on a length scale ^ r, are locally 
straight. This is the basis for the adoption of force- free initial field 
configurations in the simulations reviewed in Sec. IV. B. 

The field angular velocity, given by w = (Vi^ — V\^B^/B^)/r, is 
a field-line constant (corresponding to ip in the ideal-MHD case; see 
Sec. III. A). The ion angular velocity Vi^/r differs somewhat from uj 
but still changes only slightly along the field. Since the field lines bend 
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away from the symmetry axis, the eyhndrical radius r, and hence Vic,, 
increase along any particular field line, whereas decreases because of 
the near-Keplerian rotation law. Eventually a point is reached where 
{Vicf, — V^) changes sign. At this point, the magnetic stresses on the 
neutral gas are small and its angular velocity is almost exactly Keple- 
rian. Above this point, the field lines overtake the neutrals and transfer 
angular momentum back to the matter, and the ions start to push the 
neutrals out in both the radial and the vertical directions. It is thus 
natural to identify the base Zh of the wind with the location where 
the angular velocity of the field lines becomes equal to the Keplerian 
angular velocity. The mass outflow rate is fixed by the height Zg of the 
effective sonic point of the wind. 

As was first recognized by Wardle and Konigl (1993), one can derive 
key constraints on a viable solution by neglecting the vertical compo- 
nent of the neutral velocity, which is generally a good approximation 
throughout much of the disk column. This simplifies the problem by 
transforming the radial and azimuthal components of the neutral mo- 
mentum equation into algebraic relations. One important constraint 
that can be derived in this way expresses the intuitively obvious con- 
dition that the neutrals must be able to couple to the magnetic field 
on an orbital time scale if magnetic torques are to play a role in the 
removal of angular momentum from the disk. In the pure-ambipolar- 
diffusion regime, this condition is expressed by the requirement that 
the neutral-ion coupling time (where 7 « 3.5 x 10^^ cm^ g~^ 

s~^ is the collisional coupling coefficient) be shorter than the dynami- 
cal time. This is equivalent to requiring that the neutral-ion coupling 
parameter r] = 7Pir/Vk satisfy 

?7 > 1 . (20) 

This condition is quite general and is relevant also to disk models 
in which magnetic torques associated with small-scale magnetic fields 
transfer angular momentum radially through the disk (see chapter by 
Stone et al.; in the latter case eq. [20] identifies the linearly unstable 
regime of the magnetorotational instability, although evidently a much 
higher minimum value of r] is required for significant nonlinear growth) . 

Additional parameter constraints that can be derived by using the 
hydrostatic approximation for a wind-driving disk model in the pure- 
ambipolar-diffusion regime are given by 

(2?7)-i/2 <a<2<e?7< VK/2a, (21) 

where a = B(,/{4TTpoy/'^Cs is the ratio of the midplane (subscript 
0) Alfvcn speed to the isothermal sound speed Cg and e = — K-o/C's 
is the normalized midplane infiow speed. The first inequality corre- 
sponds to the requirement that the disk remain sub-Keplerian, the 



DISK WINDS 



21 



second to the wind-launching condition at the dink surface (subscript 
b) Br{Z\)) ^ i?o/-\/3, and the third to the requirement that the base of 
the wind he well above a density scale height in the disk. The second 
and third inequalities together imply that the vertical magnetic stress 
dominates the gravitational tidal stress in confining the disk: this is a 
generic property of this class of disk solutions that does not depend on 
the nature of the magnetic diffusivity. The last inequality expresses the 
requirement that the ambipolar diffusion heating rate at the midplane 
not exceed the rate /9o|Ko|^/2r of gravitational potential energy re- 
lease. In turn, Cs/Vk must be <C 1/(1 + e) to guarantee that the disk 
is in near-Keplerian motion and geometrically thin. One can also place 
upper limits on the density at the sonic point to ensure that the bulk of 
the disk material is hydrostatic and that M^ind does not exceed Mace- 
Similar constraints are derived by applying this analysis to the other 
diffusivity regimes identified above. The solutions that satisfy these 
constraints tend to have e ^ 1 and a }^ 1. Furthermore, the magnetic 
field in these models, whose magnitude is essentially determined from 
the condition that all the angular momentum liberated by the accreting 
matter is transported by a centrifugally driven wind, automatically lies 
in a "stability window" where it is strong enough not to be affected by 
the magnetorotational instability but not so strong as to be subject to 
the radial interchange instability (Konigl and Wardle 1996). 

As we have noted, the steady-state disk models are useful for guid- 
ing the choice of boundary conditions at the base of the outflow in 
numerical simulations of winds from quasi-stationary disks. In a com- 
plementary approach (which typically employs a highly simplified treat- 
ment of the disk outflow and therefore is not suitable for a detailed 
study of the disk wind), one can generalize the steady-state models 
and investigate the time evolution of wind-driving accretion disks (e.g. 
Lovelace et al. 1994; Konigl 1997). The time-dependent models can be 
utilized to explore the feedback effect between the magnetic flux dis- 
tribution, which affects the angular momentum transport in the disk, 
and the radial inflow induced by the magnetic removal of angular mo- 
mentum, which can modify the flux distribution through fleld-line ad- 
vection. 

C. Outflows from the Protostellar Vicinity 

As one moves to within a few AU from the protostar, the disk hydrogen 
column density can become so large (^ 192 g cm~^) that cosmic rays 
are excluded from the disk interior. (For reference, the column den- 
sity of the minimum-mass solar nebula model is 1.7 x 10'^ g cm~^ at 
1 AU). At larger columns the ionization rate is dominated by the decay 
of radioactive elements (notably ^''Al if it is present and ^°K; Stepinski 
1992), but the charge density generally becomes too low for the mag- 
netic fleld to be effectively coupled to the matter. Consequently, the 
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interior of the disk becomes inert and magnetically mediated accretion 
can only proceed through "active" surface layers that extend to a depth 
of ~ 96 g cm~^ on each side of the disk. This scenario applies both to 
a wind-driving disk (Wardle 1997) and to the case where small-scale 
magnetic fields produce an effective viscosity within the disk (Gam- 
mie 1996). Galactic cosmic rays could reach the disk along the open 
magnetic field lines that thread it, but a super- Alfvcnic outflow along 
these field lines would tend to exclude them. External ionization could, 
however, also be effected by stellar X-rays (see chapter by Glassgold 
et al.) as well as by fast particles accelerated in stellar flares. In ad- 
dition, heating by stellar irradiation could contribute to the collisional 
ionization of the surface layers (e.g., D'Alessio et al. 1998). 

The entire disk can recover an adequate coupling with the mag- 
netic field once the central temperature becomes large enough for col- 
lisonal ionization to be effective. This first occurs when the temperature 
increases above ^ 10^ K and potassium is rapidly ionized (e.g., Ume- 
bayashi and Nakano 1981). For disks in which magnetic fields dominate 
the angular momentum transport, it may be possible that Joule dis- 
sipation could maintain the requisite degree of ionization for efficient 
gas-field coupling {rj > 1; see eq. [20]). Li (1996) first explored this 
possibility for the inner regions of wind-driving protostellar disks and 
concluded that, for a disk that is in the ambipolar-diffusion regime on 
scales ^ 1 AU, a self-consistent model can only be constructed if the 
accretion rate is very high (^ 10~^ Mq yr~-^). 

A possible mechanism for recoupling the gas to the field in the 
vicinity of the YSO is the thermal ionization instability originally dis- 
cussed in the context of dwarf novae and more recently invoked as a 
possible explanation of FU Orionis outbursts (e.g., Bell & Lin 1994). 
In this picture, accretion in the innermost disk proceeds in a nonsteady 
fashion, with a "gate" at r ^ 0.25 AU opening every ~ 10^ yr or 
so after the accumulated column density has become large enough 
to trigger the instability. During the "high" phase of the instability 
(which lasts 10^ yr and is identified with an outburst) the gas is hot 
(T 10"* K) and almost completely ionized, and mass rains in at a 
rate ~ (1 — 30) x lO~^M0yr~^, whereas during the "low" phase the 
temperature and degree of ionization decline sharply, and the accre- 
tion rate drops to ^ (1 — 30) x lO~*A^0yr~^. In the context of the 
magnetized disk model, one can attribute the increase in the accretion 
rate as the gas becomes highly ionized to the reestablishment of good 
coupling between the field and the matter, which allows the field to ex- 
tract the angular momentum of the accreting gas. This could account 
both for the marked increase in Mace and for the strong disk outflow 
that accompanies it (see Sec. II. B). The magnetic recoupling idea may 
be relevant to these outbursts even if the wind is not the dominant 
angular momentum transport mechanism in the disk provided that the 
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viscosity has a magnetic origin (see chapter by Stone et al.). Wc note, 
however, that the relatively large (~ 0.1) Mwind/-^acc ratio indicated 
in some of the outburst sources is consistent with the bulk of the disk 
angular momentum being removed by the wind. 

At distances of a few stellar radii, the stellar magnetic field could 
be strong enough to drive outflows from the disk. Various scenarios 
have been considered in this connection, including individual magnetic 
loops ejecting diamagnetic blobs through a magnetic surface drag force 
(King and Regev 1994); a steady-state configuration in which mass 
is transferred to the star near the corotation radius and an outflow 
is driven along adjacent (but disconnected from the star) field lines 
(Shu et al. 1994), and time-dependent ejection associated with the 
twisting, expansion, and subsequent reconnection of field lines that 
connect the star with the disk (Hayashi et al. 1996; Goodson et al. 
1997). Numerical simulations are rapidly reaching the stage where 
they could be utilized to identify the relevant mechanisms and settle 
many of the outstanding questions. There are already indications from 
some of the existing simulations that magnetically channeled accretion 
from the disk to the star is most likely to occur when the disk carries 
a strong axial magnetic field that reconnects with the stellar field at 
an equatorial X-point; the resulting field configuration appears to be 
quasi-steady and gives rise to stiirdy ccntrifugally driven outflows along 
open fleld lines (Hirose et al. 1997; Miller and Stone 1997). 

Stellar field-driven outflows are the subject of the chapter by Shu 
et al. This class of models is based on the realization that stellar mag- 
netic field lines can be infiated and opened up through an interaction 
with a surrounding disk, and that a centrifugal wind can be driven out 
along the opened field lines. In these scenarios, the outflow typically 
originates near the inner disk radius, where the disk is truncated by 
the stellar magnetic stresses. This contrasts with the scenarios consid- 
ered in this chapter, wherein a disk-driven wind originates (and can 
contribute to the disk angular momentum transport) over a significant 
range of radii. It is worth noting in this connection that the massive 
inflows that characterize FU Orionis outbursts are expected to crush 
the respective stellar magnetospheres, so the strong outflows that are 
inferred to originate from the circumstellar disks during these outbursts 
are unlikely to be driven along stellar field lines (Hartmann and Kenyon 
1996). Coupled with the apparent inadequacy of thermal-pressure and 
radiative driving, this argument provides strong support for the rele- 
vance of disk-driven hydromagnctic winds that are not associated with 
a stellar magnetic field to the observed outfiows from FU Orionis out- 
burst sources. The ramifications of this argument become even stronger 
if a significant fraction of the mass accumulation in solar-type stars oc- 
curs through such outbursts during the early phases of their evolution 
(Hartmann 1997). To be sure, the "distributed" disk- wind models dis- 



24 



A. Konigl and R.E. Pudritz 



cussed in this chapter are meant to apply to YSOs in general and not 
just when they undergo an outburst, but only during an FU Orionis 
outburst does the disk become sufBciently luminous that an unambigu- 
ous observational signature of an extended disk wind can be obtained. 

VI. CONCLUSIONS 

Centrifugally driven winds from disks threaded by open magnetic 
field lines provide the most efficient way of tapping the gravitational 
potential energy liberated in the accretion process to power an outflow. 
The fact that such winds "automatically" carry away angular momen- 
tum and thus facilitate (and possibly even control) the accretion process 
makes them an attractive explanation for the ubiquity of jets in YSOs 
and in a variety of other accreting astronomical objects. One of the 
key findings of recent numerical simulations of MHD winds from disks 
is that such outflows are indeed easy to produce and maintain under 
a variety of surface boundary conditions. These simulations have also 
verified the ability of such outflows to self-collimate and give rise to nar- 
row jets, as well as a variety of other characteristics that are consistent 
with YSO observations. In the case of protostellar disks, the presence 
of open field lines is a natural consequence of their formation from the 
collapse of magnetically supported molecular cloud cores. Although a 
stellar magnetic field that threads the disk could in principle also play a 
similar role, this possibility is unlikely to apply to the strong outflows 
associated with FU Orionis outbursts. (More generally, it is worth 
noting that, in contrast to disk field-driven outflows, a scenario that 
invokes a stellar field may not represent a universal mechanism since 
many cosmic jet sources are associated with a black hole that does not 
provide an anchor for a central magnetic field.) Much progress has also 
been achieved in constructing global MHD disk/wind models, but the 
full elucidation of this picture and its consequences for star formation 
remains a challenge for the future. 
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